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ABSTRACT 

Four transits of the planet orbiting the star HD 209458 were observed with the STIS spectrograph on 
board HST. The wavelength domain (1180-1710A) includes Hi as well as Cl, Cil, CiV, Nv, Oi, Si, Sill, 
Sim and Siiv lines. During the transits, absorptions are detected in Hi, Oi and Cii (5±2%, 13±4.5% 
and 7.5±3.5%, respectively). No absorptions are detected for other lines. The 5% mean absorption over 
the whole Hi Lyman a line is consistent with the previous detection at higher resolution (Vidal-Madjar 
et al. 2003). The absorption depths in Oi and Cii show that oxygen and carbon are present in the 
extended upper atmosphere of HD 209458b. These species must be carried out up to the Roche lobe and 
beyond, most likely in a state of hydrodynamic escape. 

Subject headings: Star: individual (HD 209458) - Stars: planetary systems 



1. INTRODUCTION 

The extrasolar planet HD 209458b is the first one for 
which repeated transits across the stellar disk have been 
observed (~1.5% absorption; Henry et al. 2000; Charbon- 
neau et al. 2000). Together with radial velocity measure- 
ments (Mazeh et al. 2000), this has led to a determina- 
tion of the planet's radius and mass, confirming that it 
is a gas giant. During transits, the Hubble Space Tele- 
scope (HST) allowed the detection of the dense lower at- 
mosphere in the neutral sodium lines ('--^0.02% additional 
absorption, Charbonneau et al. 2002) and the extended 
upper atmosphere in Hi (^15% absorption over the stel- 
lar Lyman a emission line, Vidal-Madjar et al. 2003). This 
H I absorption extends beyond the Roche lobe, showing a 
population of escaping atoms. The escape rate is esti- 
mated to be > 10^°gs~^, consistent with theoretical eval- 
uations (Hebrard et al. 2003, Lammer et al. 2003, Lecave- 
lier des Etangs et al. 2004). The additional HST observa- 
tions which we present here allow the upper atmosphere 
of HD 209458b to be further studied. 

2. OBSERVATIONS AND DATA REDUCTION 

Four transits of HD 209458b were surveyed with the 
STIS G140L spectrograph with four HST visits on Oct. 9, 
Oct. 21, Nov. 5 and Nov. 24, 2003. For each transit, three 
exposures during three consecutive HST orbits were sched- 
uled such that at least one full exposure was obtained with 
the planet in front of the star. A previous or a following 
exposure, either entirely before the first contact or after 
the last contact, serves as a reference. 

The resulting ID and 2D spectra, output from the STIS 
pipeline, are similar to those described by Vidal-Madjar et 
al. (2003). In particular they show detectable geocoronal 
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emissions in H i and O i lines and a non-zero background 
level due to the dark current. To correct for these effects, 
we extracted the spectra from the 2D images of the de- 
tector following the procedure developed by Vidal-Madjar 
et al. (2003) and Desert et al. (2003). In each column 
of the detector we evaluate the background signal, includ- 
ing the geocorona, above and below the stellar spectrum 
and interpolate its value to be subtracted from the stel- 
lar spectrum. The errors due to the background and the 
geocorona are included in the error estimates. 

The wavelength domain of the G140L first order grating 
ranges from II80A to 1710 A. The spectra are dominated 
by various stellar emission lines (Si ill A1206A, H i A1216A, 
Nv AI239A and AI243A, Oi A1302A, Cii A1335A, Siiv 
AI394A and A1403A, Si A1474A, Sin A1527A, Civ 
A1548A and A155lA, Ci A1560A and A1657A) and the 
stellar continuum is detected down to the shortest wave- 
length (Fig. 1). The low resolution ('^2.5A) does not al- 
low the stellar emission lines to be resolved. The spectral 
information is consequently limited to the total line inten- 
sity variations. However, this allows for high sensitivity 
and large spectral coverage to search for many different 
species in the planetary upper atmosphere. 

3. ANALYSIS 

To visualize the spectral signature of the transit, we 
over-plotted off-transit and in-transit spectra (Fig. 2). A 
significant difference is seen in Hi, Oi and Cii, corre- 
sponding to an absorption during the transit. Other spec- 
tral lines with about the same intensity do not show the 
same behaviour (e.g., Sim in Fig. 2). This is the first in- 
dication that, in addition to Hi, oxygen and carbon are 
also present in the upper atmosphere of HD 209458b. 
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To quantitavely estimate the absorptions, their level of 
detection and error bars, we fitted transit profiles to the 
measured intensity of each line in each exposure as a func- 
tion of the orbital phase. The depth and width of the 
transit curve arc related to R^hs/R*, the ratio of the oc- 
culting object and stellar radii. The absorption depth due 
to the planetary disk is (i?abs/-R*)^ ~ 1.5%. The duration 
of the ingress and egress of the occupation (i.e. the width 
of the transit curve) is taken to be proportional to i?abs- 
The term i?abs/-R* is the first free parameter of the fit. 

The intensity of the emission lines shows detectable vari- 
ations from one HST visit to another. These variations 
are interpreted as stellar variations, known to occur in a 
solar type star at a few percent level on a time scale of 
several days. To extract the transit signature from these 
intrinsic variations, we consider the lines intensity as a free 
parameter for each of the 4 visits. We have thus 4 addi- 
tional free parameters for each spectral line. The intrinsic 
stellar variations arc found to be at a level of only a few 
percent between different visits. 

The stellar emissions during 3 consecutive exposures of a 
given HST visit are assumed to be constant. We consider 
unlikely that short time-scale variations could mimic a 
transit curve in H i, O i and C ii and not in the other lines, 
particularly those of the more ionized species for which 
more pronounced stellar variations should be observed. 

We measured the total line intensities after subtraction 
of the stellar continuum. The resulting 12 measurements 
have been fitted as a function of the orbital phase by an 
occultation curve with five free parameters. Results are 
shown in Fig. 3 and summarized in Table 1. The full con- 
tinuum spectrum from 1350A to 1700A shows a transit 
detection with a depth of 2.01q'7%. lines, only Hi, Oi 
and Cll show a detection at more than 2-a with transit 
depths of 5±2% in Hi, 13±4.5% in Oi and 7.5±3.5% in 
Cii (2.6, 2.6 and 2.1-(t detection, respectively). The total 
squared difference between the line intensities and the fits 
are found to have values of 8.9, 9.6 and 4.7 (for Hi, Oi 
and C II, respectively), in agreement with a distribution 
with 7 degrees of freedom. The other stellar features do 
not show detectable signatures during transit. 

A 10% depth corresponds to a radius i?abs of about 
3.6 Jupiter radii, which is also the radius of the Roche 
lobe. The measured absorption depths in O i and C ii cor- 
respond to absorbing clouds with a size of several plane- 
tary radii. The present detections show that oxygen and 
carbon are present in the extended upper atmosphere; O i 
is even present up to about the level of the Roche lobe. 

4. DISCUSSION: OBSERVATIONAL RESULTS 

4.1. Continuum 

The detection of an absorption of 2.01q"7% over the con- 
tinuum validates the adopted procedure. This absorption 
is due to the transit of the planetary disk and is known to 
be ^1.5%, well within our error bar. This shows that the 
systematic errors are significantly below this level. 

4.2. Atomic hydrogen 

Vidal-Madjar et al. (2003) detected an Hi absorp- 
tion depth of 15±4% over the Lyman-a line. The ab- 
sorption is significantly detected over only part of the 
Lyman-a line within a radial velocity range from —130 



to 100kms-i(Fig. 2 of Vidal-Madjar et al. 2003). The 
15±4% absorption in that spectral range corresponds to 
about 5.7±1.5% absorption of the total Lyman-a line in- 
tensity. The H I absorption depth measured here is consis- 
tent with the results of Vidal-Madjar et al. (2003). 

4.3. Oxygen and carbon 

The O I stellar emission triplet is a blend of the ground 

level line at 1302.2 A with the lines from excited levels 
O I* and O i** (A=1304.9 A, energy level E=158 cm'^ and 
A=1306.0A, £;=227cm-\ respectively). The Oi ground 
level line is strongly absorbed by the interstellar medium 
(Fig. 2). Therefore the ~13% absorption observed dur- 
ing the transit in the full O l triplet must be due to the 
presence of Oi* and Oi** within several planetary radii 
around the planet, up to about the radius of the Roche 
lobe. Because these excited levels of O i are likely popu- 
lated by collisions, the H i density at this altitude should 
be at least the density of a typical exobase, that is about 
lO^cm"^. Assuming solar abundances, such a density is 
needed for the Oi triplet to have enough opacity to be 
detected in a line of sight grazing the planet. In these 
conditions, only the strongest lines of the most abundant 
species in the most abundant ionization state can be de- 
tected, i.e. Hi, Oi and Cii (Ci is ionized above 11.3 eV). 

Seen at low resolution, the stellar C ii line is also a blend 
of a ground level line from C ii (A=1334.5A) with a line 
from C II* (A=1335.7 A, E=63A2 cm"!). Since the ground 
level Cii line is only partially absorbed by the interstel- 
lar medium (Fig. 2), at the present low spectral resolution, 
wc cannot a priori discriminate between absorption by C ii 
or Cii*. However, for densities above ~ 20cm~^ (Wood 
& Linsky 1997), Cii* is expected to be collisionally pop- 
ulated, we thus conclude that the detected absorption is 
due to the presence of both C ii and C II* in the extended 
upper atmosphere of HD 209458b. 

Due to the lack of spectral resolution, the real absorp- 
tion depth over a wavelength limited domain of an emission 
line (which is linked to the size of the occulting cloud) can 
be different from the measured absorption depth over the 
full stellar emission intensity (Sect. 4.2). The measured 
decrease in intensity depends on the velocity dispersion of 
the absorbing species compared to the intrinsic width of 
the stellar emission lines. HD 209458 being a solar type 
star, we can use the Sun as observed by the SUMER in- 
strument as a reference; the width of the O i and C ii stel- 
lar lines must be about ~15kms~^ for Oi and ~25kms~^ 
for Cii (Chae et al. 1998; P. Lemaire, private communi- 
cation). These values are consistent with the high resolu- 
tion echelle spectra of HD 209458 (Fig. 2). Therefore, in 
the absorbing cloud, velocity dispersions of, for example, 
25, 15 and 3kms~^ lead to corrected absorption depths 
of 13±4.5%, 13±4.5% and 65±23% for O i and 7.5±3.5%, 
12.5±6% and 62.5±30% for Cii, respectively. Note that 
we obtain similar absorption depths for Ol and Cll for 
velocity dispersions below ~15kms~^. 

5. THE ESCAPE MECHANISM 

As shown by Vidal-Madjar et al (2003), hydrogen atoms 
are escaping away from HD 209458b. This has been fur- 
ther discussed by Hebrard et al. (2003), Lammer et al. 
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(2003), Liang et al. (2003; 2004), Lecavelier des Etangs 
et al. (2004) and Yelle (2004). The present detection of 
oxygen and carbon at high altitudes provides additional 
information about the type of escape mechanism involved. 
First, oxygen is present up to about the Roche lobe where 
the Hi density must be at least 10^ cm~'^. Second, the ab- 
sorption depths of oxygen and carbon are related to their 
velocity dispersion, with ~13% depth for 15kms~^ veloc- 
ity dispersion and ~65% depth for 3kms~^. 

5.1. Jeans escape 

In the frame of Jeans escape, atoms and molecules dif- 
fuse upward from the base of the thermosphere up to the 
exobase level above which the absence of collision allows 
particles with large velocities in the tail of the Maxwellian 
distribution to escape. In the thermosphere, the scale 
height of H I is of the order of 10 000 km (0.1 planetary ra- 
dius). Under molecular diffusion conditions, the Oi scale 
height is 16 times smaller than the H i one. For an altitude 
increase of one Hi scale height, the 0/H ratio decreases by 
a factor c^*'/o. Therefore, if molecular diffusion dominates 
no 01 can be detected at high altitudes. 

However, in the solar system planetary atmospheres, 
eddy diffusion is able to transport upward species within 
the thermosphere up to some level where molecular dif- 
fusion takes over. This transport mixes the species, and 
their abundances remain unchanged upward. Assuming 
an extreme case of eddy transport up to the Roche lobe or 
to the exobase, oxygen and carbon abundances should bo 
solar at that level. However even with eddy diffusion and 
upper atmospheric temperature T^p ~10000K (Lecavelier 
des Etangs et al. 2004), thermal broadening of absorption 
lines are ~13 km for H I, and only ^3 km s~^ for O I and 
C II. This is smaller by a factor of >5 than the width of the 
stellar lines. If the absorption width was limited by the 
thermal broadening, the observed ^10% absorption would 
correspond to a geometrical absorption of more 60%. This 
would imply that oxygen and carbon extend much further 
than Hi, which is imlikcly. Eddy diffusion transport can- 
not explain the properties of the absorption in O i and C II. 

5.2. Hydrodynamic escape: "blow-off" 

To explain the present observation, we need a mecha- 
nism able to transport oxygen and carbon up to about 



the Roche lobe and to maintain their velocity dispersion 
to at least ~10kms~^. When an abundant species like 
H I escapes with a high rate, other species can be dragged 
up by the hydrodynamical flow, also called "blow-off" (see, 
e.g.. Chamberlain & Himtcn 1987). Under such conditions, 
oxygen and carbon moving out with the Hi flow should 
present a velocity dispersion at least similar to or greater 
than the sound speed, which is ~10kms~^ at 10 000 K, 
i.e. transsonic (Parkinson et al. 2004). The abundances 
of oxygen and carbon should also be preserved upward all 
along the flow. Under these conditions, the amounts of 
Oi and Cii are large enough to absorb the whole stellar 
line, allowing at least 10% absorption if they flow up to 
the Roche lobe. Although the O i and C ii absorptions are 
compatible within the error bars, the dynamical "blow- 
off" could also explain the apparent lower absorption in 
carbon as the C ii stellar emission line is broader than the 
Hi sound speed by a factor of ^2. In short, the depth 
of the O I and C ii absorptions can be interpreted through 
the absorption by oxygen and carbon flowing away in a 
hydrodynamic escape of the atmospheric hydrogen. 

6. CONCLUSION 

The still unique case of the planet HD 209458b tran- 
siting a bright star has allowed the first studies of an 
extra-solar planet atmosphere. The detection of oxygen 
(at 2.6 a) and carbon (at 2.1 cr), a signature of an atmo- 
spheric "blow-off", shows that, as already suggested by 
Hcbrard et al. (2003) and Vidal-Madjar & Lecavelier des 
Etangs (2003) hot-Jupiters can lose a significant fraction 
of their atmosphere. Ultimately, such a process could lead 
to a new type of planets with hydrogen-poor atmosphere, 
or even with no atmosphere at all (see Trilling et al. 1998). 
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Etangs 2003). Based on observations obtained with the 
NASA/ESA Hubble Space Telescope operated by the As- 
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Table 1 

Absorption depth of emission lines and continuum 



Species 


A'' 

(A) 


Absorption 

la (%) 


depth'' 
2a (%) 


Continuum 


1350-1700 






Hi 


1212-1220 






Oi, Oi*, Oi** 


1300-1310 


12.8ttj 




Cii, Cii* 


1332-1340 


7 r-l-3.6 
' -"^-3.4 





Ci 


1557- 


■1565 


4+21-1 

^•^-0.4 


<36.0 


Ci 


1654- 


■1660 


11 ^ + 11.9 
■'■■'■■"-'-11.3 


<33.6 


Civ 


1545- 


■1554 




<19.0 


Nv 


1237- 


■1246 


97 0+22.7 
^ i . o •2'j 2 


<50.0 


Si, Si*, Si** 


1471- 


•1489 


09 9+13:0 
o^./._32 2 


<58.2 


Sill, Sin* 


1525- 


1536 




<47.4 


Sim 


1204- 


■1210 


n n+2-2 

^•"-0.0 


<5.9 


Si IV 


1391- 


•1397 


fi n+6-5 
^•^-0.0 


<14.0 



^ Wavelength range for the Hnes intensity evaluation. 

^ Absorption depth given by the (i?abs/-R*)^ parameter of the fit (see 

1-a error bars and 2-a upper hmits for non-detections. 



5 




1200 1300 1400 1500 1600 1700 

Wavelength (A) 

Fig. 1. — The full HD 209458 G140L spectrum. Stellar emission lines are clearly detected together with the stellar continuum which is 
increasing towards longer wavelengths. 
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Fig. 2. — Comparison of off-transit and in-transit spectra for the Sim, Hi, Ol and Cll lines with l-cr error bars. The ofT-transit spectrum 
is the addition of the first exposure of the l*** and 3'''^ visits (thin lines); the in-transit corresponds to the addition of the fully in-transit 
exposure of each of the four HST visits (thick lines). Absorption is clearly detected in the Hi, Ol and C II lines. No signal is detected in 
other lines (e.g., Sim) with about the same amplitude. The bottom panels show the high resolution spectra obtained in 2001 with the echelle 
grating (Vidal-Madjar et al. 2003). 
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OrbiLal phase 

Fig. 3. — Plot of the lines total intensity as a function of the orbital phase. Circles, squares, triangles and diamonds are for 1^^ to 4**^ 
transits respectively. The vertical dotted lines corresponds to the position of the 1^* to 4**" contacts of the planetary disk transit. The thick 
line represents the best fit to the data (see Sect. 3). Absorptions are detected in Hi, Ol and Cll during the transits; no significant absorptions 
are detected in the other lines (i.e. Siiv). 



